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(54) A method for processing measuring values 

(57) A method for processing representations of 
measuring values, such as chromatograms, electroph- 
erograms, or spectra comprises the following steps: 

a) filtering, by a data processing means, the 
measuring values by employing a morphological or 
median filter, wherein said filter uses a structuring 
element and wherein said filtering step comprises 
the decomposition of said measuring values into a 



series of decomposed signal portions, said decom- 
posed signal portions resulting from the use of dif- 
ferent structuring elements, respectively, and 
b) analysing said decomposed signal portions with 
regard to any peaks corresponding to peaks in said 
representation of measuring values. The method 
permits precise and reproducible localization of 
peaks and determination of corresponding peak ar- 
eas, for example in a chromatogram. 
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Description 

[0001] The invention relates to a method for process- 
ing one-dimensional or two-dimensional representa- 
tions of measuring values, such as chromatograms or 
spectra or electropherograms. Such a method can be 
used, for example, to determine the peaks in a chroma- 
togram representing chromatographically separated 
sample substances. 

BACKGROUND OF THE INVENTION 

[0002] Spectra and chromatograms are often dis- 
turbed by an underlying drifting and noisy baseline. It is 
then particularly difficult to localize the peaks in the 
spectrum or chromatogram and to evaluate their area. 
This leads to inaccurate measuring results. First, errors 
in the location of a spectral or chromatographic peak 
may lead to errors in the type of substance being deter- 
mined. Second, errors in the area of the peaks lead to 
errors in the quantitative amount of the substances be- 
ing determined. 

[0003] In order to cope with the mentioned problems, 
it is known to use curve-fitting techniques, wherein ide- 
alized curves are fitted to the measured spectrum or 
chromatogram. From these fitted curves, information 
about the location and the area of the peaks can be de- 
rived. 

SUMMARY OF THE INVENTION 

[0004] in view of the prior art, it is an object of the in- 
vention to provide a method for processing measuring 
values which permits a more precise and reproducible 
localization of peaks and determination of the peak ar- 
eas. 

[0005] It is another object of the invention to provide 
a method for processing measuring values which does 
not require idealized assumptions about the shapes of 
the peaks. 

[0006] According to the present invention, these ob- 
jects are achieved by a method comprising the following 
steps: 

a) filtering, by a data processing means, the meas- 
uring values by employing a morphological or me- 
dian filter, wherein said filter uses a structuring ele- 
ment and wherein said filtering step comprises the 
decomposition of said measuring values into a se- 
ries of decomposed signal portions, said decom- 
posed signal portions resulting Irom the use of dif- 
ferent structuring elements, respectively, and 

b) analysing said decomposed signal portions with 
regard to any peaks corresponding to peaks in said 
representation of measuring values. 

[0007] According to the invention, it has been found, 
surprisingly, that morphological or median filters can be 
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employed advantageously for analyzing one-dimen- 
sional or two-dimensional representations of measuring 
values, in particular chromatograms or spectra. In ac- 
cordance with an embodiment of the invention, the ortg- 

5 jnal measuring values, lor example the chromatogram, 
are decomposed into various signal portions, wherein 
different decomposed signal portions correspond to dif- 
ferent structuring elements used for the decomposition. 
From the various decomposed signal portions, the pre- 

10 cise locations and the areas of the peaks can be deter- 
mined. The invention can also be used for analyzing 
electropherograms, i.e. : representations of measuring 
values obtained by an electrophoretic separation proc- 
ess, such as capillary electrophoresis. 

is [0008] In a preferred embodiment for analyzing chro- 
matograms, the structuring element is an elliptical struc- 
turing element and the width of the structuring element 
is changed successively during various iterations of the 
filtering process, such as to produce different decom- 

20 posed signal portions. 

[0009] The method of the invention has several ad- 
vantages. A first advantage is that only few parameters 
are needed. According to an embodiment of the inven- 
tion, a very good approximation of the ideal signal can 

25 be achieved with only three estimated parameters: a) 
minimal height, b) minimal width, and c) maximal width. 
[0010] A second advantage of the method of the in- 
vention is its robustness. The method can treat nearly 
any kind of signal on the basis of only very few estimated 

30 parameters. For example, the method of the invention 
yields accurate and linear results for signals with nearly 
any kind of baseline and signal-to-noise ratio. In this 
connection, the term "accurate" is to be understood in 
terms of relative error, and the term "linear" is to be un- 

35 derstood with regard to variation of the height and area 
of the peaks. According to the invention, the baseline 
drift can be eliminated completely, noise is drastically 
reduced and the evaluation of the peaks is satisfactory. 
[0011] A third advantage of the method of the inven- 

40 tion is that it does not deform the peaks. It substantially 
preserves the foots, shoulders or any other features of 
the peaks without any prior assumptions about shapes 
being required. Unlike the classical applications of 
mathematical morphology, the present invention does 

4S not just use mathematical morphology to recognize pat- 
terns or to detect objects, but it actually extracts peaks 
and thus enables a further, more detailed analysis, such 
as separation of overlapping peaks, determination of 
peak areas, and other analyses. 

so [0012] Furthermore, the method of the invention inte- 
grates with a filtering part a part for precisely localizing 
and evaluating the areas of the peaks. Thus, it is not 
necessary to use additional specialized software and 
the results become more accurate. The invention thus 

55 provides for automatic recognition of the peaks in a 
chromatogram or spectrum while respecting the areas 
of the peaks. 

[0013] The morphological filters (i.e., non-linear fil- 
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ters) used in the invention are representing information 
in terms of shapes and not in terms of frequency as most 
of their standard linear counterparts. This leads to a bet- 
ter preservation of sharp edges and renders the method 
of the invention especially appropriate for peak recog- 
nition. Thus it offers in many cases (such as edge de- 
tection, noise removal) a better alternative to the clas- 
sical frequency based approach. Moreover, such filters 
are easily adjustable to the data dimensions: one-di- 
mensional signals, images or higher dimensional sig- 
nals can be treated equally well. The mentioned mor- 
phological filters offer a great versatility. Due to the pos- 
sibility to choose the reference form, i.e. the structuring 
element, and due to the virtually unlimited possibilities 
of combining basic operators, one can customize the fil- 
ters to treat a wide range of problems. 
[001 4] An important aspect of the present invention is 
the discovery that special types of mathematical filters 
can be used in chromatography or spectrophotometry 
for signal analysis, yielding unexpected advantages. 
Some types of mathematical filters are described, in a 
different context than the present invention, in the fol- 
lowing publications: 

J. Serra, "Image Analysis and Mathematical Morpholo- 
gy", London (1981); D. Wilson et aL, "Morphological 
Restoration of Atomic Force Microscopy I mages Lang- 
muir 1995, 11, 265-272. 

Another prior art method for processing measuring val- 
ues is known from: G. Pearson, "A General Baseline- 
Recognition and Baseline-Flattening Algorithm", Jour- 
nal of Magnetic Resonance, 27, 265-272 (1997). 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] In the following, embodiments of the invention 
will be explained in detail with reference to the drawings. 

Figure 1 illustrates some operations used in mor- 
phological filters . 

Figure 2 shows the effect of morphological opera- 
tions on a signal. 

Figure 3 illustrates results achieved with a method 
of the invention in analyzing chromato- 
grams. 
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togram typically consists of a collection of data points 
representing the variation of a detector signal measured 
at the outlet of the chromatographic separation column 
as a function of time. The data points are processed by 
5 a digital data processing means which may be integrat- 
ed with the liquid chromatograph or which may be a sep- 
arate unit connected to the chromatograph. The meas- 
ured and the processed chromatograms as well as peak 
locations and peak areas resulting from the processing 
10 can be displayed on a display screen and/or printed out 
by request of a user. When using the invention in con- 
nection with a spectrophotometer, the measuring values 
comprise a collection of data points representing inten- 
sities at different wavelengths. The spectrophotometer 
* s can be used, for example, to perform a spectral analysis 
of the substances eluting from the separation column of 
a liquid chromatograph or of a capillary electrophoresis 
device. 

[001 7] The method of the invention makes use of mor- 
phological filters. Some morphological filters have been 
developed in morphological mathematics which is a the- 
ory for analyzing data with regard to form and structure. 
Morphological filters are non-linear filters representing 
information in terms of shape and not in terms of fre- 
quency as most of the standard linear filters. With mor- 
phological filters, the analysis is performed with the help 
of an adapted element, the "structuring element". An in- 
troduction to morphological mathematics can be found 
in the above mentioned textbook by J. Serra. In the fol- 
lowing, some terms of morphological mathematics 
which will be used in the subsequent description will be 
explained. 

[0018] Basic operators in morphological mathematics 
are dilation ©, erosion ©, closing and opening. For the 
explanation of these operators, reference is now made 
to Figure 1A-C. In Figure 1A-C, reference sign M de- 
notes the signal which is to be filtered, and L denotes 
the structuring element with a reference point marked 
by a dot. The umbra of a signal is defined as all the val- 
ues ■under" the signal surface, i.e. the grey part in Figure 
1A. 

[0019] The dilation M©L of the signal M by the ele- 
ment L is defined by the following formula: 

M©L = max. , (M(x-j)+L(j)) (x-j) eM 
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Figure 4 illustrates the localization of peaks in a 
chromatogram from the pre-processed sig- 
nal, so 

Figure 5 illustrates the principle of a sieve filter. 

DETAILED DESCRIPTION OF THE INVENTION 

55 



[0016] The method of the invention deals with the 
processing of measuring values, such as a chromato- 
gram produced by a liquid chromatograph. The chroma- 



[0020] The dilation operation can be understood 
graphically (Figure 1 C) as placing the reference point of 
the element L on all points within the umbra of M. L thus 
describes an output umbra which defines the dilated sig- 
nal M®L 

[0021] The erosion of the signal M by the structuring 
element L , i.e. the operation M0L is described by the 
following formula: 

M©L = min jeL (M(x + j)-L(j)) (x+j) eM 
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[0022] The erosion operation can also be understood 
as moving vertically , for each point, the structuring el- 
ement L such that it remains entirely in the umbra of M. 
This is illustrated in Figure 1B. 

[0023] The aforementioned two basic operators can 
be combined to obtain the closing operator (•) and the 
opening operator <°) as illustrated in Figures 2A and 2B. 
M denotes the original signal and L the structuring ele- 
ment L. The dotted curve in Figure 2A is the result of the 
opening operation, and the dotted curve in Figure 2B is 
the result of the closing operation. The opening opera- 
tion as illustrated in Figure 2A consists in first eroding 
and then dilating the signal M with the same structuring 
element L: 

M (o) L = M0L0L 

[0024] The closing operation as illustrated in Figure 
2B consist in first dilating and then eroding the signal M 
with the same structuring element L: 

M (•)!_= M©L0L 

[0025] The opening operation can also be illustrated 
by the following analogy. The operation corresponds to 
a process wherein one takes an object whose top sur- 
face is described by the structuring element, presses the 
object upwards against the underside surface of the sig- 
nal and moves it along the signal. Finally, the curve de- 
scribed by the highest positions of the object corre- 
sponds to the opened signal. Thus, positive features of 
the signal that are smaller than the structuring element 
are eliminated. Figure 2A shows how the structuring el- 
ement L is pressed, according to the mentioned analo- 
gy, from below in the triangular and rectangular ele- 
ments of the signal. 

In case of a spherical structuring element, the opening 
operation is called "rolling ball algorithm". 
[0026] In the following, the term "sieve" will be ex- 
plained. A morphological filter can roughly be consid- 
ered as a method for extracting all the features from a 
signal which are, first, either positive or negative and, 
second, wider than the structuring element. To get a full 
lengthy description of a signal it is easy to imagine that 
one may decompose it on various levels, each corre- 
sponding to a fixed size of the features. This is called a 
sieve. It is obtained by the following method: 

1 . Filter the signal S n .-, with the structuring element 

2. Subtract the filtered signal S n from the original 
signal to gain the residues R n . 

3. Increase the width of the structuring element (re- 
specting its form) from one (n = n + 1 ) and iterate 
the process. 



The sieve based on the filter \j/ n (n being the width of the 
structuring element) is defined as follows; 

S m( A ) = V m <- V 2 (Vi(A)...) 

Where A is the original signal. The residues are calcu- 
lated according to the formula: 

R m = S m ., (A)-S ra (A) 

[0027] The collection of residues (Rj) (i = 1 ,2 : ...,n) plus 
possibly the last filtered signal Sn + 1 then represent the 
signal with its features classified amongst the residues 
according to their width. 

[0028] Figure 5 illustrates a sieve filter with an exam- 
ple. The original signal to be filtered is labeled with ref- 
erence numeral 50. Block 51 depicts a sieve decompo- 
sition of the signal 50 using a morphological filter. The 
various levels 1 /£ , 1 , 2, 3, 4, 5 in block 51 are generated 
by varying the width of the structuring element 52. The 
widths of the structuring element 52 are schematically 
illustrated at the right of each level. The width of the 
structuring element is successively increased from the 
first to the last level. As can be seen, the original signal 
50 is thus decomposed into various levels, from which 
the peaks in the signal are clearer apparent and from 
which quantitative analyses of the peaks can be made. 
[0029] A useful property of sieve filters is the relatively 
small influence of noise on the final result. Typically, 
high-frequency (or short-term noise) is for the greatest 
part concentrated on the very first levels of the sieve, 
corresponding to the small structuring elements. There- 
fore, it can be easily removed from the final result, more- 
over, it does not unduly influence the other levels of the 
sieve. 

[0030] The sieve filters are also able to compactly de- 
scribe and efficiently classify the different features con- 
tained in the original signal. In contrast to matched filters 
or wavelet analysis, these filters are binding the features 
with edges. That means that peaks, baseline or any no- 
ticeable feature would be found in a few number of levels 
and always at the same place, limited by the edges of 
the selected feature. 

[0031 ] The method of the invention is particularly use- 
ful for analyzing chromatograms, as will be described in 
the following. The data set representing the chromato- 
gram is first submitted to an opening morphological filter 
with an adapted structuring element. The structuring el- 
ement is adapted with regard to the minimal height of 
the expected peak. First, the structuring element is cho- 
sen so as to be very small in terms of length. The filter 
lets pass only the elements contained in the data set 
which are of larger width than the structuring element. 
The filtered curve is subtracted from the original data 
set, thus giving the first residue. The process is then it- 
erated with the elements growing in length. As a result, 
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the residues are describing on the first levels the high- 
frequency noise, which can also be described as small 
non-smooth parts of the signal. Then, the chromato- 
graphic peaks appear, and on levels corresponding to 
structuring elements wider than the maximal peak width, 
drift and baseline of the chromatogram can be found. 
The useful signal is easily gained back by adding the 
residues of the middle levels. In this connection, "middle 
levels" corresponds to a size of the structuring element 
between the minimal and maximum width of the peak. 
[0032] In the filtered signal, the locations of the peaks 
are then roughly estimated with the help of the minimal 
height parameter This parameter can typically be cho- 
sen by a user such that it is greater than the noise level, 
but smaller than the smallest peak. Any value in this 
range leads to good results, but it is preferred to keep 
some distance from the limits ol the range. 
[0033] The underlying peak area is then evaluated by 
precisely locating the peaks in the residues of the sieve 
and adding the involved surface. The peak localization 
in this embodiment has been performed in two steps, i. 
e., by an initial coarse localization and a subsequent fine 
localization, in order to achieve localization faster and 
more reliably. 

[0034] Figure 3 illustrates in more detail the results 
which can be achieved with the invention in connection 
with chromatograms. The horizontal axis is the time ax- 
is, the indicated numbers are data points. The vertical 
axis represents the intensity. 

The original curve, i.e., the chromatogram 35, is shown 
at the top of Figure 

3. The ideal signal is denoted with reference numeral 
36, and the signal resulting from the filtering is denoted 
with reference numeral 37. Together with the filtered sig- 
nal 37, a curve 38 is shown , which indicates the devia- 
tion between the filtered signal 37 and the ideal signal 
36. 

[0035] In the embodiment illustrated in Figure 3 : an 
elliptical structuring element 31 is used. On the various 
levels of the sieve, the height 32 of the structuring ele- 
ment is kept constant while the "radius" or width 33 is 
decreasing throughout the levels of the sieve. An ellip- 
tical structuring element has certain advantages for an- 
alyzing chromatograms because there is typically a 
great difference between vertical and horizontal resolu- 
tion in chromatograms. The operator used in this em- 
bodiment is a basic opening operator as defined above. 
[0036] Processing the levels of the sieve: 
The range of levels (i.e., the range of widths of the struc- 
turing elements corresponding to each level) used to re- 
build the signal is critical to eliminate undesired effects 
without distorting the peaks. The minimal peak width 
and the maximal peak width provides the range of levels 
that are to be selected. The minimal peak width is de- 
termined by the resolution of the chromatographic ap- 
paratus and the maximal peak width is a parameter to 
be chosen as being smaller than the wander and base- 
line variations. The maximal peak width is typically cho- 
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sen by the user. The other levels typically contain only 
short term or long term noise effects and can be 
dumped. Within that range, levels are to be correctly 
spaced, that means that the increase of the structuring 

5 element is dimensioned such that the levels will contain 
a satisfying image of the peaks. In a practical example, 
for peaks varying between 24 and 250 data points in 
length, the "radius" of the elliptical structuring element 
is made to vary between 10 and 300 data points in steps 

10 of 10 for a height corresponding approximately to 10 
times the minimal peak height. 

[0037] A further improvement can be made by intro- 
ducing a threshold in the levels. This removes minor 
contributions which are most of the time due to noisy 
is effects. The criterion can be the median value of the con- 
sidered level, with a threshold value proportional to it. 
[0038] When side effects are present, a further im- 
provement can be achieved by widening the signal be- 
yond its start and end points by extending its end values 
20 for the maximal size of the structuring element. This is 
illustrated in Figure 4, wherein the ranges 45 and 46 rep- 
resent the extended end portions. 
[0039] With the sieve method of the invention, the er- 
ror is minimal and a significant improvement is made in 
25 terms of reproducibility and linearity as compared to 
known filters. Furthermore, the sieve filter is particularly 
resistant to noise. Other useful features are the flatness 
and the absence of horizontal offset. 
[0040] Through the sieve process, access is gained 
30 to a lot more information than actually used just for elim- 
inating unwanted noises. In particular, by analyzing the 
levels, a precise location of the peaks can be extracted. 
Based on this knowledge, a simple addition of the peaks' 
values will give a good approximation of the area under 
35 the peaks. That means that a complete analysis of the 
signal is performed, completing the first filtering and 
avoiding the need of further software. 
[0041] The following algorithm can be used for inte- 
grating the peak areas. For explaining the algorithm, ref- 
40 erence will be made to Figure 4, wherein the upper part 
shows a filtered signal and the lower part shows the ex- 
tracted peaks. 

1 . Points that belong to the peaks according to the 
45 criterion "value is greater than half the minimal peak 

height" are marked in the reconstructed signal (Fig- 
ure 4, threshold line 41). 

2. For each marked point, the adjoining areas of 
non-zero terms are searched in the various levels. 

so The values found are added to form a new signal. 

3. The non-zero parts of the newly gained signal are 
separated, each corresponding to a peak, hereby 
giving precisely the locations of the peaks. This is 
illustrated in the lower part of Figure 4 by the dots, 

55 such as dots 42 and 43. The underlying area is de- 
termined by adding the value of each point of the 
concerned peak. In the example shown in the lower 
part of Figure 4, the first peak extend between 
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points 42 and 43. In terms of data points (horizontal 
axis), this corresponds to a range between data 
points between number 472 and 534. The values at 7. 
each of these data points are added to obtain the 
integrated peak area. s 

[0042] Various modifications to the above described 
embodiment are possible. For example, different mor- 
phological filters having other structuring elements 
(such as a segment or a Gausssian element) could be to 
used. Also, more elaborated operators can be used. In- 
stead of a morphological filter, a related but non-mor- 
phological filter, such as a median filter could be used. 
[0043] Modifications are also possible with regard to 
choice and use of the levels. For example, the space is 
and range of the levels may vary. Furthermore, other 
methods to identify the peaks, for example a gradient 
method can be used. 



is an operator performing an opening operation 

A method as in any of the preceding claims, wherein 
different decomposed signal portions (52a, 52b, ...) 
are generated by varying the width of the structuring 
element in a predetermined way 



Claims 
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1. A method for processing one-dimensional or two- 
dimensional data sets representing or derived from 
measuring values, such as chromatograms, elec- 25 
tropherograms, or spectra, 
characterized by the following steps: 



a) filtering, by a data processing means, said 
data sets (50) by employing a morphological or 
median filter, wherein said filter uses a structur- 
ing element and wherein said filtering step com- 
prises the decomposition of said data sets into 
a series of decomposed signal portions, said 
decomposed signal portions (52a : 52b, ...) re- 
sulting from the use of different structuring ele- 
ments, respectively, and 

b) analyzing said decomposed signal portions 
with regard to any peaks corresponding to 
peaks in said data sets representing or derived 
from said measuring values. 



30 



35 



40 



3. 



4. 



5. 



A method as in claim 1 , comprising a step of deter- 
mining the positions of said peaks in said decom- 
posed signal portions. 45 

A method as in claim 1 or 2, additionally comprising 
a step of determining the area of said peaks in said 
decomposed signal portions. 



A method as in any of the preceding claims, wherein 
said filter is a morphological filter comprising a 
structuring element and an operator. 

A method as in claim 4, wherein said structuring el- 
ement is an elliptical structuring element. 



so 
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6. A method as in claim 4 or 5, wherein said operator 
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